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Abstract 

We present calculations on the magnetocapacitance of the two- 
dimensional electron gas in a lateral semiconductor superlattice under 
two-dimensional weak periodic potential modulation in the presence 
of a perpendicular magnetic field. Adopting a Gaussian broadening 
of magnetic- field-dependent width in the density of states, we present 
explicit and simple expressions for the magnetocapacitance, valid for 
the relevant weak magnetic fields and modulation strengths. As the 
modulation strength in both directions increase, beats of the magneto- 
capacitance oscillations are observed, in the low magnetic field range 
(Weiss-oscillations regime), which are absent in the one-dimensional 
weak modulation case. 



KEYWORDS: lateral surface superlattices; density of states; magnetocapaci- 
tance. 
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1 Introduction 



A lateral surface superlattice (LSSL) combines a system of nano-fabricated 
gates with a MODFET, where a periodic pattern is imposed onto 2DEG at 
the semiconductor heterostructurep]. This pattern can be a one- or two- 
dimensional periodic potential modulation of different strength allowing for 
a variety of artificial periodic structures. The periodic potential is applied 
through an array of metal gates whose bias can be varied[2]. The modu- 
lation which is purely induced by the modulated gate potential, leads pre- 
dominantly to a sinusoidal potential shape and its strength can be tuned by 
a nano-patterned top gate as well as by a back gate electrode [3]. Magne- 
totransport measurements in LSSL devices have revealed novel oscillations 
in the conductivity, the Weiss oscillations [1] which are superimposed on top 
of the well-known Shubnikov-de Haas oscillations. Weiss oscillations were 
first discovered in LSSL devices with holographically defined one-dimensional 
(ID) superlattices[S], and were explained^ E] in terms of a commensu- 
rability criterion between the cyclotron orbit diameter at the Fermi level 
2R C = 2^2im e l 2 (where n e is the electron density and I = \J h/eB the 
magnetic length) and the superlattice period a: 

2R c =(\ + <p)a A = 0,1,2,... (1) 

where is a phase factor. 

Besides Weiss oscillations, novel magneto- resistance oscillations with 1 /B 
periodicity have been observed in short-period superlattices when one flux 
quantum h/e passes through an integral number of lattice unit cells[8] 126] . 
A large fraction of these studies was performed on square pH] or rectangular 
lattices [TTl[T2"] . while only a few experiments on hexagonal lattices have been 
published P31 E]. For the majority of the effects, the lattice type is, in 
principle, irrelevant. Nevertheless, Altshuler-Aronov-Spivak oscillations have 
been observed around zero magnetic field in hexagonal lattices, while Aronov- 
Bohm oscillations can be detected in larger magnetic fields [Sj. 

The capacitance of a mesoscopic structure is of electrochemical nature. It 
depends in an explicit manner via the thermodynamical density of states on 
the electronic properties of the structure. The investigation of the magnetic- 
field dependence of the electrochemical capacitance plays a significant role 
in several experiments, especially those using capacitance spectroscopy. Ca- 
pacitive studies of the density of states of a two-dimensional electron gas are 
possible because the thermodynamical density of states appears as a series 
capacitance to the geometrical capacitance [IS]. The behavior of the elec- 
trochemical capacitance can be also be examined from a dynamical point of 
view[16j. At low frequencies, the frequency dependent capacitance can be 
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expanded in frequency and at the linear order it is a series combination of 
the static capacitance and a charge relaxation resistance. This approach has 
been performed to the analysis of mesoscopic capacitors in Refs fTTl [T8] . 

The oscillatory characteristics caused by commensurability between the 
cyclotron diameter 2R C and the modulation period a in a LSSL, reveal 
themselves in the electrochemical magnetocapacitance measurements. A 
pronounced modulation of both the minima and maxima of the capaci- 
tance oscillations has been observed in LSSLs with ID weak periodic po- 
tential and explained as a consequence of the oscillatory bandwidth of the 
modulation-broadened Landau levelsjl]. Although in recent years special at- 
tention has been paid to investigate the magnetic miniband structure as well 
as the magneto-transport properties in a LSSL under two-dimensional weak 
modulation [SI EH], magnetocapacitance calculations pertinent to this case are 
rather limited [T9]. 

In this paper we present calculations on the density of states and the 
magnetocapacitance of the 2DEG in a LSSL under two-dimensional weak 
periodic potential modulation in the presence of a perpendicular magnetic 
field. Adopting a Gaussian broadening of magnetic-field-dependent width, 
we present explicit and simple expressions for the DOS, valid for the relevant 
weak magnetic fields and modulation strengths. We investigate the magne- 
tocapacitance, as a function of the magnetic field, for various strengths of 
the periodic potential modulation. In the low-field range (Weiss-oscillations 
regime), a beating structure in the oscillating magnetocapacitance is ob- 
served. 

The paper is organized as follows: In Section 2 we briefly present the 
magnetic spectrum of the 2DEG in a LSSL under two-dimensional weak 
periodic potential modulation. In Section 3, the DOS is calculated analyti- 
cally and numerically as a function of energy and magnetic field. In section 4 
we present numerical results concerning the magnetocapacitance oscillations. 
We conclude with a summarizing Section 5. 

2 Energy spectrum under weak 2D-periodic 
modulation potential 

To describe the electrons in the conduction band of a lateral superlattice 
at the AlGaAs-GaAs interface in a constant perpendicular magnetic field 
B = Bz, we employ a model of strictly 2DEG with a 2D periodic potential 
modulation V(r), r = (x,y). Following the common practice, we adopt an 
effective mass approximation to take into account the effect of the crystalline 
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atomic structure over the charge carriers. The system is described by the 
well-known single-particle Hamiltonian of Bloch electrons in a homogeneous 
magnetic field: 

H = ^( P + eA) 2 + V(r) (2) 

where we neglect the Zeeman splitting and spin-orbit interactions which are 
very small in GaAs systems. The effective mass m* for electrons in AlGaAs- 
GaAs is taken m* = 0.067 m e , where m e is the free-electron mass. 

We consider a two-dimensional periodic modulation with rectangular sym- 
metry, i.e: 

2it 2tt 
V(x, y) = V x cos ( — x) + V y cos( — y) (3) 

d x dy 

where V x , V y are the modulation strengths and a x , a y the periods in the 
corresponding directions. We assume that the modulation is weak compared 
to the cyclotron energy hw c so that mixing of different Landau levels can 
be neglected. The homogeneous magnetic field is represented by the vector 
potential A = (-By/2, -Bx/2, 0). 

Applying a direct perturbation theory [26] or a projection operator approach [20] 
to unmodulated system, the energy spectrum is obtained as a function of the 
crystal momentum (k x , k y ) and the number of flux quanta through a lattice- 
unit cell $/$o, where $o = h/e is the magnetic flux quantum. 

E(k x , k y ) = (n + - )hu c + V x F n (u x ) cos (— u x k y ) + V y F n (u y ) cos (— u y k x ) 

I 71 7T 

(4) 

where F n (u) = exp (-u/2)L n (u), u x = (ir$ /$)a y /a x , u y = (-K® /<f>)a x /a y 
and L n (u) are the Laguerre polynomials. The perturbation approach has 
been shown to be an extremely good approximation to the direct diagonal- 
ization of the Hamiltonian, for the parameter values of interest (V x ,V y << 
Ep, where Ep is the Fermi energy) and for magnetic fields 0.1 T < B < 
1 T. Thus, the unperturbed Landau levels broaden into bands, with width 
W n (B) = 2{y x \F n {u x )\ + Vy\F n (u y )\), that oscillates with magnetic field B. 

In the low-magnetic-field range it is a good approximation to take the 
large n-limit of the Laguerre polynomials [2*Tj that appear in W n (B). For a 
square lattice with the width of the bandwidth is then obtained 

as 

(2t\R c TT 



WJB) = HW±£ 



7T 




COS 



- 
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(5) 



where R c = lyf2np + 1 is the cyclotron radius with n F = E F /hu c the Landau 
index at the Fermi level. From equation (JS]) we obtain that the bandwidth 
reaches zero (flat band) when 2R c /a = (A+l/4) and maximum when 2R c /a = 
(A + 3/4) with A = 0,1,2,... . 
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3 Density of States 



The DOS per unit area of the 2DEG, in the absence of disorder, can be 
expressed as a series of 5-functions: 

°w = ^E^ £ - ( 6 ) 

Tikx ky 

where g s = 2 is the electron-spin degeneracy. 

In the case of ID electrostatic modulation, the 5-functions in Eq. §§§ 
are broadened into bands with Van Hove singularities at the edges of each 
Landau band. On the other hand, if 2D modulation is imposed, the DOS does 
not exhibit Van Hove singularities reflecting the 2D nature of the electron 
motion in the Landau bands with finite mean velocities |20j. 

In practical two-dimensional electron systems, there is always some broad- 
ening present due to scattering centers. We assume a Gaussian broadening 
of zero shift and a width T which, according to short-range scattering theory 
of Aoki and Ando[22], is proportional to \f~B. In this case Eq. ([6]) takes the 
form 



2tt 3 / 2 



98 E [ dt * [ ^T^f ex P ( " ^ - E ntxty f/2T 2 ^ 



(7) 



We expand the integrand of Eq. (j^J) in powers of e n)ti = ViF n (ui)cos(ti) , 
i = x,y, assuming weak modulation strengths. The odd powers of e njii vanish 
after integrating over t x and t y . Then, to second order in the modulation 
potential, the even powers can be easily integrated and finally we obtain 

^7 = 7^f^ e 1 + — fs (8) 



where A 2 n = \V x F n (u x )\ 2 + \V y F n (u y )\ 2 and D = m*/irh 2 is the DOS of a free 
2DEG at B = 0. Near the center of the nth Landau level, the expression is 
valid only when \ViF n {ui)\ 2 <T,i—x,y. 

Using the Poisson's summation formula [2"Tj. we carry out the summation 
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over n and the DOS becomes 

2 



D(E) 



D 



o 



1 -^{-^jK x exp(-(n x r/2E) 2 ) sm(2K x ) + cos(2n x )/K x 



1 fV \ r 
+2^ (-l) n exp ( - (mrr//k,; c ) 2 ) cos {2imE/hw c ) 



(9) 



2vm 2 / V, \ 2 



fkOr 



COS [K x r 



4' 



2vm 2 / K \ 2 



tlUJr 



COS r 



4' 



where 



^3 



2E 2nl 



2E 2nl 



(10) 



Figure 1(a) shows the DOS as a function of energy. We have adopted 
a field-dependent width T = 0.35y/B meV which is consistent with most 
of the earlier numerical studies in the low magnetic field range [23] . The 
other parameters used in the numerical calculation are V x = V y = 0.3 meV, 
B = 0.64 T, a x = a y = 365 nm. The solid curve shows the result from the 
analytic expression ([9]) and the doted one is obtained numerically from Eq. 
(J7J). As it can be seen, the agreement between the curves is good for most 
of the energies. In Figure 1(b) we show the DOS at the Fermi energy as 
a function of the magnetic field. The Fermi energy is evaluated for given 
electron density n e inserting Eq. in the relation 



dEf(E - E F )D(E) 



'111 



where fix) = [1 + e x ^ KT ]~ 1 is the Fermi-Dirac distribution function. For the 
parameters listed above, the correction to the Fermi energy due to weak two- 
dimensional modulation, is found to be of order 10 -3 over the range B < 0.7 T 
and hence, for n e = 3.2 x 10 11 cm~ 2 , a fixed value of Ep = 11.43 meV has 
been adopted. Good agreement is found between the numerical and the ap- 
proximate curves except at certain magnetic fields for which the amplitude 
of the DOS-oscillations is slightly underestimated. The DOS exhibits a peak 
at each Landau-band center. We observe that the weak 2D-periodic poten- 
tial produces clear modulation of the envelope of the DOS-oscillations and 
a beating structure appears in the range 0.15 T < B < 0.45 T (Weiss - 
oscillations regime). The above beating structure is not observed in the ID 
weak modulation case. 
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Figure 1: Calculated DOS of a 2DEG in a periodic 2D electrostatic mod- 
ulation ( function of energy at B = 0.64 T and (b) as a function 
of magnetic field B. The solid curve corresponds to exact numerical re- 
sults and the dotted to approximate from Eq.Q. The parameters used are 
V x = V y = 0.3 meV, a x = a y = 365 nm, E F = 11.43 meV, T = 0.35^ meV. 

4 Beats in the magnetocapacitance oscilla- 
tions 

The capacitance of a system consisting of a metal gate- insulator-semiconductor 
sandwich (e.g. gated AlGaAs/GaAs heterostructure), depends not only on 
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the thickness of the insulator but also on the DOS at the semiconductor side 
and on the material's parameters. If the two depletion layers interpenetrate 
each other, the gate voltage Vg, is connected to the electron density by [24"] 

V G = -n e + ^ + K (12) 

€i e 

where L is the thickness of the AlGaAs-layer, £j is the dielectric constant 
of the layer and K is a constant that takes into account fixed charges in 
the AlGaAs. Differentiating Eq.( fT2l with respect to n e , one obtain the total 
inverse magneto-capacitance [25] per unit area C(B) at a given temperature 
T: 



1 1 1 



D T (0) 



C(B) C e 2 D T (0) L D T (B) 
where Co denotes the total capacitance at zero magnetic field, and 



(13) 



Expression (I13p is valid when a change in the gate voltage affect only 
the charge in the 2DEG and the gate, leaving intact the charge of impu- 
rities in the heterostucture. This condition is fulfilled at low temperatures 
in a LSSL based on AlGaAs/GaAs heterostructure. Following a dynami- 
cal approach [TO], an analogue expression has been recently derived by Wang 
et.a/.[2S] for the frequency dependent electrochemical capacitance in order 
to study high-frequency inductive corrections in a quantum capacitor similar 
to the experiment of Gabelli et.al.[27\ 

Figure 2 shows our numerical results for the capacitance between the gate 
and the 2DEG as a function of magnetic field at temperature T = 1.5 K. We 
consider the case of square-symmetric modulation with a x = a y = 365 nm 
and V x — V y — 0.1, 0.2, 0.3, 0.4 meV. The zero-filed capacitance has been 
taken from experiment [29J equal to Co = 1168 pF/ mm 2 . As the modulation 
strength in both directions increases, modulated magnetocapacitance oscil- 
lations with nodes (a beating pattern) are observed in the low magnetic filed 
range 0.15 T < B < 0.45 T. The phases of the oscillations changes at these 
nodes and the number as well as the amplitude of beatings increases as the 
modulation strength increases. We should note that the potential strengths, 
are kept in the weak modulation regime (V x « Ep, V y « Ep), so that Lan- 
dau Level mixing is prevented. Concerning larger magnetic fields B > 0.5 
T (Shubnikov-de Hass oscillations range), the envelope of magnetocapaci- 
tance maxima increases monotonically with increasing magnetic field, while 
the the envelope of the magnetocapacitance minima decreases monotonically 
with increasing magnetic field. 
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Figure 2: Calculated magnetocapacitance versus magnetic field for potential 
strengths V x = V y = 0.1, 0.2, 0.3, 0.4 meV, scattering broadening T = 
0.35v^B meV, and periods a x = a y = 365 nm. 



The origin of the above beating structure is the oscillatory behavior of 
the bi-directionally modulated bandwidth. This behavior is clearly shown 
in Figure 3 where the contribution 5C(B) = C(B) — Co is plotted together 
with the bandwidth oscillations in the magnetic- field range of interest. The 
beatings of the magnetocapacitance oscillations, at a certain magnetic field 
B, correspond to adjacent minima or maxima of the modulated bandwidth 
W n (B) and the amplitudes of the oscillations under each envelope are closely 
related to the value of W n (B). In other words, the beatings occur nearly 
in the middle between adjacent 'flat-band' conditions. Especially, for the 
beating structure of Figure 3, we found that the cental peaks of the beatings 
appear at the magnetic fields B = 0.157 T, 0.185 T, 0.227 T, 0.292 T. These 
values coincide with the magnetic fields given by the following commensura- 
bility relation derived from Eq.([T]) 



B, 



ea(A + 0) 



(15) 



S 




Figure 3: Modulated bandwidth W n (B) and beating patterns of magneto- 
capacitance oscillations for the square LSSL with a x = a y = 365 nm and 
V x = V y = 0.4 meV. 

with = ±1/4 and A = 2, 3. 



5 Concluding Remarks 

We have studied the magnetocapacitance oscillations of a 2DEG in a LSSL 
under 2D weak periodic modulations at low temperatures (T ~ 1 K). Adopt- 
ing a Gaussian broadening of magnetic-field-dependent width, an explicit an- 
alytic expression for the oscillatory behavior of the DOS has been obtained. 
An overall agreement between the exact and approximate results is found 
for the relevant weak magnetic fields (B < 1 T) and modulation strengths 
(V x « Ep, V y « Ep). The calculated magnetocapacitance has been shown 
to have a rich oscillatory structure in this regime. As the modulation strength 
in both directions increases, a beating pattern is observed in the low mag- 
netic field range 0.15 T < B < 0.45 T due to the oscillatory behavior of 
the bi-directionally modulated Landau-level bandwidth. It is our aim to 
extend our calculations for the case of a square LSSL with non-symmetric 
2D modulation as well as the case of a LSSL with different periods in both 
directions. We are not aware of any directly relevant experimental work. 
We hope though that the results described above will motivate experiments 
in which the low-field magnetocapacitance could be measured in a weakly 
bi-directionally modulated 2DEG of a LSSL. 
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